Interplay between lattice-scale physics and the quantum Hall effect in graphene 
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Graphene's honeycomb lattice structure underlies much of the remarkable physics inherent in this material, 
most strikingly through the formation of two "flavors" of Dirac cones for each spin. In the quantum Hall regime, 
the resulting flavor degree of freedom leads to an interesting problem when a Landau level is partially occupied. 
Namely, while Zeeman splitting clearly favors polarizing spins along the field, precisely how the states for each 
flavor are occupied can become quite delicate. Here we focus on clean graphene sheets in the regime of quantum 
Hall ferromagnetism, and discuss how subtler lattice-scale physics, arising either from interactions or disorder, 
resolves this ambiguity to measurable consequence. Interestingly, such lattice-scale physics favors microscopic 
symmetry-breaking order coexisting with the usual liquid-like quantum Hall physics emerging on long length 
scales. The current experimental situation is briefly reviewed in light of our discussion. 

PACS numbers: 73.43.-f, 71.10.-w, 71.10.Pm 



I. INTRODUCTION 

Whereas electrons confined to 2D in conventional GaAs 
heterostructures are well-represented as residing in a con- 
tinuous space, the behavior of electrons in graphene is in- 
extricably tied to their underlying honeycomb lattice struc- 
ture. This simple fact is central to much of the exceptional 
physics in the graphene quantum Hall effect. Most dramat- 
ically, in zero field the honeycomb band structure supports 
two Dirac cones centered at momenta ±Q = ±(47r/3, 0), 
so that graphene realizes a "relativistic" dispersion (over a 
significant energy range). In the presence of a perpendicu- 
lar magnetic field, the Dirac spectrum develops into Landau 
levels that are approximately four-fold degenerate, one fac- 
tor of two coming from spin and the other from the two "fla- 
vors" of Dirac cones. Taking into account this degeneracy 
as well as a Berry phase acquired from the cyclotron motion 
leads to the prediction of quantized Hall plateau s with 
a xy = ve? jh = 4(j + l/2)e 2 /h, whose observation was one 
of the early achievements in the study of graphene^. 

The Landau level structure as well as the formation of 
these plateaus can be readily captured by a continuum Dirac 
equation with SU(4) symmetry (reflecting invariance under 
spin and flavor rotations), with no further considerations of 
graphene's lattice structure. In a sense, the problem then be- 
comes rather similar to GaAs 2D electron systems, with a 
modified spectrum. But is this the full story? Simple sym- 
metry considerations suggest that it is not. The approximate 
four-fold Landau level degeneracy noted above is clearly bro- 
ken down by Zeeman coupling, but this still leaves an SU(2) 
flavor symmetry in the continuum theory. This residual sym- 
metry can certainly be spontaneously broken, as occurs due to 
Coulomb exchange in quantum Hall ferromagnetism wherein 
plateaus occur at all integer filling factors (see the reviews 
of Refs. 0Ht). However, it is important to observe that there 
is no microscopic reason for its existence in the first place. 
This suggests that something can, and will, explicitly lift this 
degeneracy further. Such effects were absent in initial exper- 
iments, due presumably to their being overshadowed by dis- 
order. Given the observation of additional v = 0, ±1, and 



±4 plateaus at higher magnetic fields 910 as well as the in- 
evitable future enhancement of sample quality, the question 
of precisely how the flavor degeneracy is lifted in graphene 
quantum Hall states is not only interesting, but also relevant. 
Such issues are important to resolve for a number of reasons. 
For instance, what is the precise nature of these new integer 
quantum Hall states? Looking forward, what type of states 
should be anticipated in the fractional regime? How do the 
low-lying excitations behave? What, if any, symmetries are 
broken in such states? And how can such properties be de- 
tected? 



The central goal of this short review is to demonstrate 
that subtler lattice-scale physics provides such a degeneracy- 
lifting mechanism. There are a number of possible sources 
for this physics, which have been briefly summarized in Ref. 
|7|. Here we will expand on some of these ideas, concentrat- 
ing on clean graphene systems in the regime of quantum Hall 
ferromagnetism, which we believe to be of current experimen- 
tal relevance. For the most part, we will examine a continuum 
Dirac theory derived from a simple lattice Hamiltonian that in- 
corporates the electron kinetic energy and Coulomb repulsion, 
and demonstrate that such a minimal model indeed propagates 
important flavor-symmetry-breaking interactions into the con- 
tinuum. Interestingly, this model predicts that quantum Hall 
states at filling factors v = ±1 are not liquid-like at all on 
the lattice scale, but rather exhibit microscopic charge den- 
sity wave order. Similarly, a broken-symmetry pattern with a 
tripled unit cell is favored at other odd-integer v. (See Fig. [2]) 
We then briefly summarize a recent numerical study of such an 
interacting lattice model 11 , which confirms these conclusions 
in the clean limit and extends the analysis by incorporating 
disorder. An alternate degeneracy-lifting mechanisms origi- 
nating from disorder-driven variations in the electron hopping 
strengths as proposed by Abanin et alM is also reviewed. 
Finally, we conclude with a brief assessment of the current 
experimental situation regarding these subtle yet interesting 
lattice-scale effects. 
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II. CONTINUUM THEORY FROM AN INTERACTING 
LATTICE MODEL 

A minimal lattice model for graphene that can be expected 
to realize integer and fractional quantum Hall states when a 
strong magnetic field is turned on should consist of two parts: 
electron kinetic energy plus Coulomb repulsion. It is reason- 
able, then, to start by understanding the physics contained in 
such a Hamiltonian, which we write using second quantiza- 
tion as 



H 

Hke 



H, 



Coul — 



Hke + H Coui , 

-* J2 J2 [4xw+h.c], 

(xx') Q = T,1 
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In these expressions, c^ x is the electron creation operator for 



spin a, n ax = c' a 
ber operator, n x = 



is the corresponding electron num- 



ij x + ri| x , and S(x) 



/3C/3x 



is the usual spin operator with er a vector of Pauli matrices. 
The first term in HqouI contains the on-site repulsion energy, 
written in a manifestly SU(2)-invariant form, while the second 
represents the long-range part of the Coulomb repulsion, with 
V(x.) the Coulomb potential. Normal ordering of the opera- 
tors in Eq. (01 is understood. 

The kinetic energy term in the Hamiltonian can be readily 
diagonalized to reveal two Dirac points at momenta ±Q = 
±(47r/3,0), which he at the Fermi level when there is one 
electron per site. To describe excitations near the two nodes, 
it is highly convenient to pass to a continuum formulation by 
expanding the momentum-space fermion operators in terms 
of two flavors of continuum Dirac fermion fields (denoted R 
and L) as follows, 



~ 7^aiia(q) 



(4) 
(5) 



where 7 = ^/^/(S 1 / 4 !) (I is the linear system size). Here and 
below we reserve indices a, (3 for spin, A, B for flavor, and 
a, b for the honeycomb sublattice. It is also useful to employ 
Pauli matrices <7 3 a/3 , t ab , and rf ab that contract with the spin, 
flavor, and sublattice indices, respectively. We will use the 
convention that suppressed indices on the fields are implicitly 

summed {i.e., Ipty = J2 a Aa ^aAa^Aa)- 

With these identifications, it is straightforward to recast the 
full interacting Hamiltonian into a continuum Dirac theory. 
Upon turning on the external magnetic field, one can essen- 
tially guess the leading kinetic energy and interaction terms, 
namely 

Ho = -ihv J <Px^[r) x D x +r) y D v ]i) 

+ iy"d 2 xd 2 x , p tot (x)U(x-x')ptot(x'), (6) 



with v ~ 10 6 m/s the Fermi velocity, Dj = dj — i(e/K)Aj, 
and ptot = the total density. For concreteness we take 
the magnetic field B = V x A in the +z direction, normal 
to the graphene sheet. These contributions clearly exhibit an 
SU(4) symmetry, being invariant under arbitrary spin/flavor 
rotations. A more careful derivation, however, reveals a num- 
ber of anisotropy terms, which break this SU(4) symmetry 
down to U(l) sp in x [U(l) x Z 2 ]fl aV or- Writing the full Hamil- 
tonian as H = Ho + Hi, these additional interactions can be 
written as follows, 

Ul = j ~ S^bB ■ Stat + ^"o[Pt 2 ot + Pstag 

l(S 2 R1 + S 2 L2 + 6S m -S L2 + (1^2))] 



^2 Vl ( r Wag(x + r)p stag (x) 
r 

u 2 [jIj+ + 



(7) 



where Stot = \i> '(Tij) is the uniform spin density, /O s tag = 
ip<T z r] z ip represents the staggered electron density between 
sublattices 1 and 2 of the honeycomb (see Fig. Syia = 

\ijr Aa crijjA a , and J + = ip' m tpL2 and J_ = ip^tpLi represent 
components of the density oscillating at wave vectors ±2Q. 
The sum in Eq. © is over Bravais lattice vectors. Once again, 
normal ordering is understood in the interactions above. 

While the presence of these short-range anisotropy terms is 
not a priori obvious (apart from the Zeeman coupling), each 
has a simple physical interpretation. The second term in Eq. 
O merely encodes the on-site U repulsion; displaying the lat- 
tice constant ag explicitly, we have 



u 



= V3alU/A. 



(8) 



The third, whose coupling is given by 



Vi(t) 



°[V(r + a /V3y)-(l-5 TtO )V(r)}, (9) 



represents an inter-sublattice repulsion that reflects the smaller 
Coulomb energy cost for electrons residing on the same sub- 
lattice versus opposite sublattices. While this interaction is 
certainly short-ranged, the above non-local form must be re- 
tained when dealing with filling factors v = ±1 to have an 
effect. At other integer filling factors, however, a purely local 
form suffices, and one can replace 

wi(r)p sta g(x + r)p s tag(x) -> uip 2 tag , (10) 



It! 



1 J e 2 



(11) 



where e is the (unscreened) dielectric constant. The final term 
represents the intra-sublattice repulsion between density com- 
ponents oscillating at ±2Q, with 



u 2 « 



(12) 
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We note that the u-i interaction was also noticed by Goerbig 
et al., 13 though at odd-integer filling factors we find that the 
vi term provides the leading flavor-symmetry-breaking inter- 
action. 

Before turning to the implications of these interactions 
for quantum Hall states, it is useful to note the hierarchy 
of energy scales in the problem. The kinetic energy term 
in the SU(4)-invariant Ho gives rise to four-fold degener- 
ate Landau levels 1,14 indexed by an integer n, with energies 
E n = sign(n)y/2ehv 2 B\n\. The spacing between the n = 
and n = ±1 levels is the largest energy scale at roughly 
400 v / £ [T] K, where B[T] is the magnetic field evaluated 
in Teslas. The characteristic Coulomb energy scales with the 
field in the same way: 



£c 



47T6 



RPA^B 



100-v/-B[T] K, 



(13) 



where £b is the magnetic length and erpa is the screened 
dielectric constant computed within the random phase 
approximation 15 , which yields erpa ~ 5eo- The energies as- 
sociated with the symmetry-breaking terms in Hi are much 
smaller at laboratory fields. The Zeeman energy for instance 
is gusB ~ -B[T] K. Furthermore, as emphasized in Ref . 1 1 3l 
the energies for the short-range interactions in Hi are down 
by factors of aojls compared with £c- The associated scale 
for these terms is 

^short-range = f ^ ~ B[T] K, (14) 

the first factor being the characteristic lattice-scale Coulomb 
energy and the second reflecting the average number of elec- 
trons in the highest occupied Landau level per unit cell. As 
we will discuss below, an exception occurs at filling factors 
v = ±1, where the scale for flavor symmetry breaking is 
down by an additional factor of ao/is- Given the clear sepa- 
ration in energy scales above, it is reasonable to first consider 
the physics contained in Ho, and then take into account the ad- 
ditional symmetry-breaking terms in Hi. This is the strategy 
we will follow below. 



in. IMPLICATIONS FOR QUANTUM HALL STATES 

Let us now very briefly review some aspects of integer 
quantum Hall physics expected from the SU(4) invariant 
Hamiltonian Ho- Apart from the integer quantum Hall states 
occurring at filling factors v = A(j + 1/2), where the highest 
occupied Landau level is completely full, it is well-established 
that quantum Hall ferromagnetism emerges at other integer 
filling factors due to Coulomb exchang o 7 i 8 i 16 . Here, in the 
absence of anisotropy terms the SU(4) symmetry enjoyed by 
Ho is broken spontaneously, giving way to gapless Goldstone 
modes analogous to spin wavesi^ (the number depends on 
v) and SU(4) skyrmions 18 ' 19 that provide the lowest-energy 
charge excitations in the n = 0, . . . , ±3 levels^. 

The physics of quantum Hall ferromagnetism in graphene 
is substantially modified by the anisotropy terms in Hi. Only 



a subset of the SU(4) symmetry will be spontaneously bro- 
ken, the Goldstone modes will generally become gapped, and 
the skyrmion character will change as well. Here we will fo- 
cus for the most part on elucidating the qualitative features 
of these effects, obtained within a Hartree-Fock framework; 
quantitative aspects have been discussed elsewher o 17 ' 20 i 21 . We 
will specialize to the n = and n = 1 Landau levels (n = — 1 
is related by particle-hole symmetry). To this end, we first 
note a simple yet crucial feature of the Landau level wave 
functions that will prove illuminating. Namely, the n = 
single-particle wave functions for flavor L reside entirely on 
honeycomb sublattice 1, while flavor R states reside entirely 
on sublattice 2. In contrast, the probability weight for n ^ 
wave functions is evenly distributed between both sublattices. 
The effect of Hi in these two cases differs completely as a 
result of this distinction. 

Consider the n = Landau level first. Employing the 
standard Landau level projection, the interactions simplify 
greatly here due to the character of the wave functions noted 
above, with the U2 term dropping out entirely. At filling factor 
v = — 1, corresponding to a quarter-filled level, Zeeman cou- 
pling clearly favors a spin-polarized state; the on-site U inter- 
action can then play no role due to Pauli exclusion. The fate of 
the flavor degree of freedom is thus set by the sublattice repul- 
sion in Hi- It is useful to note that a general flavor-polarized 
state can be expressed as 



n 

m=0 



i<t>J 
2" "\L,m 



vac) 



(15) 



where cl A m adds a spin-up, flavor- A particle with angular 
momentum m into the n = Landau level and |vac) corre- 
sponds to the empty n = level. The angles 9, <fi specify the 
polarization direction in "flavor space", shown schematically 
in Fig. Q] In particular, "easy-axis" flavor polarization is fa- 
vored here, with 9 = or it, resulting in a state with all the 
electrons occupying the n = level spontaneously choosing 
to reside on one honeycomb lattice as shown in Fig.|2ja). The 
physics behind this microscopic charge density wave (CDW) 
is simple: the electrons can stay farther apart from one another 
by remaining on one sublattice. 

It is important to note that this effect is rather weak. First, 
the gap Af for exciting "flavor waves" (which weakly modu- 
late the system away from the perfect CDW, much like a spin 
wave) is quite small: 



4 x 10 _3 (B[T]) 3/2 K, 



(16) 



which is down by an extra small factor of ao / is compared to 
^short-range in Eq. ( TBI . Despite appearances, a local p 2 tag 
interaction is insensitive to 9 and <f) above, for the same reason 
that (S z ) 2 is trivial for a spin-1/2 moment. Rather, it is the 
non-local part of the sublattice repulsion in Hi that provides 
the degeneracy lifting, resulting in the additional small factor 
noted above. The spin-wave gap, by contrast, is much larger 
at A s = gusB w B[T] K; hybrid "spin-flavor waves" cost 
a similar energy. The smallness of the flavor anisotropy im- 
plies that "flavor skyrmions" will provide the lowest-energy 
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FIG. 1: Order-parameter space for flavor-polarized states. Interac- 
tions favor "easy-axis" polarization at v = ±1, but "easy-plane" 
polarization at v = 3, 5. The corresponding lattice-scale order is 
shown schematically in Fig.[2] Disorder of the type discussed in Sec. 
V, however, favors easy-plane polarization at v — ±1. 



(a) (b) 




FIG. 2: Lattice-scale symmetry-breaking patterns favored by inter- 
actions at (a) v = ±1 and (b) v = 3, 5. On the left, the electron den- 
sity is higher on one of the two sublattices. On the right translation 
symmetry is broken, along with, in general, rotation and reflection 
symmetries due to the formation of lattice-scale currents circulating 
around hexagons. For example, one such ordering involves currents 
circulating uniformly around the shaded hexagons. 

charge excitations, but also suggests that other mechanisms 
not included here may influence the flavor degree of freedom, 
as we will discuss below. Second, only the relatively small 
number of electrons in the n = level participate in the CDW 
(around 1.4 x 10~ 5 £>[T] electrons per hexagon), though in 
principle its presence may be detected with STM or NMR 
measurements^. Analogous results hold for the particle-hole 
related v = 1 filling factor. 

The ground state at v = depends on the competition be- 
tween Zeeman coupling and the on-site U versus the sublattice 
repulsion, which is much more effective here since even a lo- 
cal /9g tag interaction behaves nontrivially on polarized states. 
If the former dominate, then the ground state will clearly be 
spin-polarized, with electrons equally occupying both sublat- 
tices. This transpires provided gpsB + 2uoPo > 4:UiPo, 
where po = l/(2nl 2 B ) is the density of a quarter-filled Lan- 
dau level. For clean graphene sheets on a SiC<2 substrate, we 
estimate that this inequality is satisfied if U > 4 eV. Other- 
wise the sublattice repulsion will dominate, and a spin-singlet 



CDW emerges. We emphasize that here the CDW order can 
be much more robust than at v = ±1, the flavor-wave gap 
now being set by Auipa — (gpsB + 2uopo) which should 
generically be of order £ short-range- Which scenario prevails 
is unclear because of uncertainties in the magnitude of U 22 
though a spin-polarized state can always be achieved by tilting 
the field to enhance the effective ^-factor. In either case, the 
charge gap is expected to be set by mixed spin/flavor-textured 
skyrmions. For instance, in the spin-polarized state skyrmions 
will have lower energy if they degrade the spin polarization 
while simultaneously restoring CDW order. The dependence 
of the activation energy on an in-plane field can allow one to 
distinguish between these ground states. An increase is ex- 
pected in the spin-polarized state since the larger effective g- 
factor reduces the optimal skyrmion size and hence raises their 
energy, while a decrease is expected for the spin-singlet CDW 
since here skyrmions induce a revival of spin polarization. 

Next we discuss the n = 1 Landau level, beginning with 
filling factor v — 3. The physics is similar for v = 5 and 
will not be discussed separately. Since the wave functions 
now live equally on both sublattices, one must wrestle with 
all terms in 7i\. Satisfying the Zeeman coupling requires 
full spin polarization, so again the on-site U repulsion drops 
out. A general flavor-polarized state can be written as in Eq. 
( fT31 ), with the creation operators now filling n = 1 Landau 
levels. Both the u\ and U2 terms lift the ordering degener- 
acy here, though the effects are subtle since the form of the 
wave functions dictates that (p s tag) = (J±) = for any 
9, <fi. Rather, the symmetry-breaking comes from exchange 
contributions. In particular, ui favors "easy-plane" polariza- 
tion, i.e., 9 = 7r/2 and arbitrary <j), while U2 favors easy-axis 
polarization. A quantitative estimate suggests that the former 
easy-plane state, which is characterized by an order param- 
eter (i/j'(t x + iT y )tp) 0, has lower energy. This state 
is ordered at wave vectors ±Q, and exhibits the symmetry- 
breaking pattern shown in Fig. |2|b). In general, rotation and 
reflection symmetries are broken as well, due to the forma- 
tion of lattice-scale currents circulating around honeycomb 
plaquettes. The specific pattern of currents is determined by 
(f>, which within our theory is selected spontaneously. Con- 
sequently, one would expect a finite-temperature Kosterlitz- 
Thouless transition, which has yet to be definitively observed 
in a quantum Hall system. While graphene may eventually 
provide a clean setting for realizing this physics, we note that 
(presumably) higher-order lattice effects will energetically fa- 
vor specific values of cf>. 

At filling factor v = 4, the ground state depends on the 
strength of Zeeman coupling and the on-site U relative to the 
ui and u-2 interactions. The former terms favor a fully spin- 
polarized state, while the latter favor a spin-singlet state with 
easy-plane flavor polarization. Variational energetics suggest 
that for reasonable values of U the ground state is spin polar- 
ized, and hence does not exhibit lattice-scale order. 

Similar physics can be expected at higher- \n\ Landau lev- 
els, although this problem has not yet been examined quanti- 
tatively. With regard to fractional quantum Hall states, which 
have been studied theoreticall y 3, 17 i 23 i 24 i 25 i 26 but not yet ob- 
served, the anisotropy terms considered here should be kept 
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in mind, particularly in situations where a number of states 
are energetically competitive. 



IV. NUMERICS 

The picture presented above for the nature of quantum Hall 
states at odd-integer filling factors is supported by recent nu- 
merics by Sheng et alM- These authors carried out an exact di- 
agonalization study of the interacting Hamiltonian in Eq. {T), 
incorporating the magnetic field directly on the lattice (i.e., 
sending t — ► te lA ™' , where A xx > is the lattice vector poten- 
tial, and including the Zeeman energy) and additionally allow- 
ing for random on-site chemical potential disorder. Assuming 
the Zeeman splitting is sufficiently large so that spin-up and 
spin-down Landau levels are well-separated, the many-body 
wave functions were obtained exactly upon projecting onto 
the highest occupied Landau level. Up to 24 electrons in the 
projected level were considered, in systems with linear dimen- 
sions up to 200 x 200 and periodic boundary conditions. 

In the clean limit appropriate to the discussion above, spin- 
and flavor-polarized ground states are found numerically at 
filling factors v = — 1 and v = 3. In particular, a clear 
easy-axis anisotropy is observed at v = —1, corresponding 
to the lattice-scale CDW order of Fig. Ufa). The anisotropy 
energy characterizing the energy difference between the easy- 
axis and easy-plane polarized states (see Fig. [TJ is indeed 
found to be quite small, on the order of £c(o,o/^b) 2 as ex- 
pected from our Hartree-Fock analysis. Again, we note that 
this is down by a factor of ao/lg compared to £ s hort-range 
since a local p 2 tag interaction does not split the ordering de- 
generacy. Unfortunately, the smallness of the anisotropy en- 
ergy did not permit numerical resolution of the correspond- 
ingly small gap A f for flavor-waves out of the easy-axis state. 
At v = 3, easy-plane flavor polarization emerges. A strong 
finite-size effect supports the lattice-scale structure predicted 
for this state [see Fig.[2fb)]. Specifically, the easy-plane state 
occurs only if both linear dimensions of the system are divisi- 
ble by 3; otherwise easy-axis polarization occurs, since the or- 
dering at wave vectors ±Q would then be "frustrated" at the 
boundaries. Finite-size scaling suggests that the easy-plane 
state wins in the thermodynamic limit, since the ground-state 
energy per particle exhibits an upturn when the easy-axis state 
occurs. Thus, numerics are entirely consistent with our ana- 
lytic predictions for the clean limit, and, importantly, can pro- 
vide valuable information about the stability of quantum Hall 
ferromagnetism when disorder is included, as we will mention 
briefly at the end of this review. 



V. DISORDER EFFECTS AT v = ±1 

We now return to filling factors v = ±1, where the flavor 
symmetry breaking encoded in the interacting lattice model 
we have been focusing on was found to be rather weak. 
Abanin et al. 12 have proposed an interesting alternative ef- 
fect which may provide the dominant flavor symmetry break- 
ing mechanism — disorder. Real graphene sheets inevitably 



exhibit some degree of structural imperfection, which may 
arise for instance from the presence of a substrate. Atomic 
force microscopy on graphene samples examined in Ref. 
revealed surface ripples with a typical lateral length £ of a 
few tens of nanometers (though ensuing fabrication advances 
produced higher-mobility samples with ripples below the res- 
olution of their microscope 2 ^). Such distortions give rise to 
spatially varying electron hopping strengths, which appear as 
an effective inhomogeneous magnetic field in the continuum, 
directed oppositely for the two flavors. In our notation this 
corresponds to sending 

Dj -► Dj - i(e/h) aj (x)T z (17) 

in Eq. ©, where 6h(x) = V x a is the effective field, whose 
magnitude was estimated to be around 0. 1 to 1 T in Ref. |27L 
The vector potential a, being proportional to the local hopping 
strength deviations, was assumed to have white noise correla- 
tions with a correlation length £ in Ref. 

One possible consequence of this random field is that the 
system may break up into domains with local easy-axis fla- 
vor polarization [i.e., local CDW order) specified by the di- 
rection of Sh, at the cost of domain wall energy. This is the 
Larkin-Imry-Ma stat o 28,29 . Abanin et al., however, argue that 
an easy-plane flavor-polarized state is energetically more fa- 
vorable, at least for weak disorder. The physics behind this 
easy-plane selection can be understood by analogy to a ferro- 
magnet, where spins preferentially align perpendicular to an 
applied field so they they can effectively gain energy by relax- 
ing toward it. Just as we discussed at v = 3 and 5, the easy- 
plane state here will exhibit a Kosterlitz-Thouless transition; 
the transition temperature was estimated to be a few Kelvin in 
Ref. [HI The easy-plane anisotropy arising from disorder was 
estimated to be 

\BJ \£bJ 10y/B[T\ 

where Sh[T] is the typical effective field strength in Teslas and 
£[nm] is the disorder correlation length in nanometers. With 
5h in the range of 0. 1 to 1 T and £ ~ 30 nm, this may indeed 
significantly exceed the intrinsic anisotropy A f arising from 
interactions. For systems with less structural imperfection 27 , 
however, the competition between disorder and interactions 
may be quite delicate. 

VI. CONCLUDING REMARKS 

We conclude by discussing the current experimental situ- 
ation regarding the physics of quantum Hall ferromagnetism 
in graphene. Although v = ±4 plateaus have been observed 
at high magnetic fields (B > 20 T), 910 quantum Hall ferro- 
magnetism is apparently not yet realized in the n = 1 Landau 
level. Plateaus at v = 3,5 have yet to be resolved, and the 
activation energy at v = 4 was found to be dominated by 
Zeeman splitting of the spin-up and spin-down states, rather 
than Coulomb energy as would be expected for a clean sys- 
tem. The situation is entirely different for the n = Landau 
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level, where v — and ±1 quantum Hall states have all been 
observed*^. Furthermore, the activation energy measured at 
v = 1 was several times larger than that at v = ±4 and found 
to scale with \f~B, suggesting Coulomb exchange as the origin 
of the gap 10 . Thus current experimental graphene samples ap- 
pear sufficiently clean to enable observation of quantum Hall 
ferromagnetism in the n = level, but not at other levels. 
It is worth emphasizing that the greater stability of ferromag- 
netism to disorder in the n — level is fully consistent with 
theory. Nomura and MacDonald derived a Stoner criterion 
for the onset of interaction-driven quantum Hall states, which 
predicts a critical mobility for spontaneous symmetry break- 
ing that is several times larger in the n — 1 level than in the 
n = level 16 . More recently, the numerical study by Sheng 
et al. u finds a critical on-site disorder strength for observing 
a plateau at v — 3 to be roughly one third that at v = 1. 

Resolving the more detailed symmetry-breaking order in- 
herent in the likely exchange-driven v = and ±1 quantum 
Hall states remains an open experimental challenge. Ulti- 
mately, local probes such as STM and NMR measurements 



would be ideal for detecting possible microscopic patterns 
such as those in Fig. [2] though such experiments may be diffi- 
cult since the order would involve only a small fraction of the 
total number of electrons. Such symmetry breaking may also 
be revealed by possible concomitant lattice distortions. On a 
cruder level, the qualitative trend in the activation energy as 
a function of field might provide some clues as to the nature 
and origin of these states. At v = ±1 for instance, the in- 
trinsic anisotropy due to interactions scales as B 3 / 2 , while the 
anisotropy due to disorder predicted by Abanin et al. scales 
as B^ 1 / 2 . Consequently, the energy for skyrmions relative to 
those in the SU(4)-symmetric case should increase with field 
in the former case, but decrease in the latter, which may be 
revealed through measurements of the activation energy. 
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